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We study the tunneling density of states (DoS) in the superconducting systems driven by Zeeman 
splitting Ez into the paramagnetic phase. We show that, even though the BCS gap disappears, 
superconducting fluctuations cause a strong DoS singularity in the vicinity of energies —E* for 
electrons polarized along the magnetic field and E* for the opposite polarization. The position of 



this singularity E* = | \Ez + \J E 2 Z — A 2 J (where A is BCS gap at Ez = 0) is universal. We found 

analytically the shape of the DoS for different dimensionality of the system. For ultra-small grains 
the singularity has the shape of the hard gap, while in higher dimensions it appears as a significant 
though finite dip. The spin-orbit scattering, and the orbital magnetic field suppress the singularity. 
Our results are qualitatively consistent with recent experiments in superconducting films. 

PACS numbers: 74.40. +k, 74.50. +r, 73.40Gk, 73.50.-h 



I. INTRODUCTION 

It is well known that the magnetic field, H , suppresses superconductivity since it lifts time reversal symmetry (see 
e.g. Ref. ^ for general introduction). In the absence of the spin orbit coupling, this effect can be separated into two 
mechanisms: (i) effect of the magnetic field on the orbital motion associated with Aharonov-Bohm phase, and (ii) 
Zeeman splitting of the states with the same spatial wave functions but opposite spin directions. 

In the bulk systems, the suppression of the superconductivity is typically associated with the first mechanism. 
Indeed, the estimate for the critical field H C2 in this case is 

H C2 e - 0o, (i.i) 

where <j>o = t£ is the superconducting flux quantum and 

(1.2) 

is the coherence length for the dirty superconductors, A is the BCS gap, and D is the diffusion coefficient. On the 
other hand, the magnetic field necessary to affect the superconductivity by virtue of the spin mechanism is given by 



gLHBH spin ~ A, (1.3) 



where is the Lande <?-factor, and /i# = is the Bohr magneton. Comparing Eqs. (1.1) and ([Lg|), one finds that 
H spin is far in excess of H C2 : 

%^e F r»l, (1.4) 

-nc 2 



where ef is the Fermi energy and r is the elastic momentum relaxation time. Condition (1.4) means that in the bulk 
system, the orbital effect of the magnetic field is always dominant . 

Situation may change in the restricted geometries. Consider, e.g., the superconducting film of the thickness a <C £, 
placed in the magnetic field parallel to the plane of the film. The Cooper pair in this case is restricted in the transverse 
direction by the film thickn ess a . As a result, the geometrical area swept by this pair can be estimated as at; rather 



than as £ 2 . Therefore, Eq. (IT) should be changed to 



ijj>-0o => Hl~H C2 [i). (1.5) 



On t he o ther hand, Zeeman splitting Ez — gLl^BH is not affected by geometrical restriction. Accordingly, instead of 



Eq. (1.4) the ratio of the two scales of magnetic field is given by 



1 



(1.6) 



Thus, for sufficiently thin films, a <C -^p, the spin effects become dominant. One can easily check that the same 
estimate (1.6) holds for other restricted systems, i.e., superconducting grains or wires. In these cases, a is the size of 
the grain or the diameter of the wire respectively. Quite generally, a is determined by the minimal size of the sample 
in the plane perpendicular to the magnetic field. In this paper we consider such restricted geometries, and unless the 
opposite is stated, neglect the orbital effects. P 

The transition from superconductor to paramagnet is of the first ordera: superconducting state is the only stable 
state at Ez < A; while at Ez > 2A normal state is the only stable state. Both phases are locally stable in the 
interval of magnetic fields where A < Ez < 2A. The normal state becomes lowest in energy and thus globally stable 
at Ez > v2A. From now on, we will assume that this condition is fulfilled. 

One of the most fundamental manifestations of the superconductivity is the gap in the tunneling density of states 
(DoS) around the zero biasou. One can expect that after the paramagnetic transition not only BCS order parameter 
vanishes but also the energy dependence of tunneling DoS becomes similar to those in superconductors above critical 
temperature T c . (The latter dependence is discussed in the review article Ref. g) 

In this paper, we demonstrate that, on the contrary, there are clear observable superconducting effects in the normal 
state even far from the transition region. We will show that at the transition point there appears a dip in the DoS 
(schematic evaluation of the DoS with the magnetic field is shown in Fig. 0). 
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FIG. 1. Evolution of the tunneling DoS with the Zeeman 
splitting Ez for (a) the superconducting state, Ez < %/2A, 
see e.g. Ref. |^, and (b) for the paramagnetic state E z > \/2A. 
The usual zero bias anomaly in the paramagnetic state (b) 
is not shown for simplicity. Shape of the singularity at E* 
corresponds to the O-dimensional case. 



The shape and the width of this dip depend on dimensionality of the system. However, its position is remarkably 
universal: 



E* = \{E z + ^El-&), 



(1.7) 
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for OD (grain), ID (strip), and 2D (film) cases. Some of conclusions have been already briefly reported by two of usB. 
Here we present detailed derivations of results of Ref. [| and consider the relevant perturbations (spin-orbital coupling, 
orbital magnetic field, finite temperature, and energy relaxation) of the new tunneling anomaly. 

The remainder of the paper is organized as follows. Section || presents the parametrically exact solution for the 
simplest but instructive geometry of zero-dimensional systems (ultra-small superconducting grains). Section III deals 
with more involved problem of the tunneling anomaly in the superconducting films and wires. Both sections required 
application of the diagrammatic technique on the level of at least Ref. 0. For the benefit of the readers interested 
in physical interpretation rather than in rigorous derivations, we present in Sec. IV the qualitative derivation which 
grasps all the essential physics involved, even though fails to give the completely quantitative description. Section |v| 
analyzes how the tunneling anomaly is affected by spin-orbital coupling, orbital magnetic field, finite temperature, 
and energy relaxation. We discuss the recent experimenia on the Zeeman splitting of the tunneling anomaly in Al 
films in Sec. VI. Our findings are summarized in Conclusion. 



II. ZERO DIMENSIONAL SYSTEMS 



Let us consider an isolated disordered superconducting grain which is small so that the Zeeman splitting dominates 
over the orbital magnetic field effect (see, e.g, Refs. p|pO| for recent experiments on such grains). We assume that the 
size of the grain exceeds electronic mean free path I, and, at the same time, it is much smaller than the superconducting 
coherence length £. We also assume that kpl ^> f . This results in the large dimensionless conductance of the grain 
g (g ~ kpla). Finally, we assume that the grain is already driven into the paramagnetic state by the Zeeman splitting. 
Our goal is to find effects of the superconducting fluctuations on the DoS of the system. 

The Hamiltonian H of the system consists of noninteracting part Hq and interacting one Hint ■ Using the basis of 
the exact eigenstates of Ho labeled by integers i and j, one can write the Hamiltonian as 

H = ^E la a\ a a. la - \Sj2 a l] a li a jl a j1- (2- 1 ) 

icr i.j 

Here operator a\ a (ai a ) creates(annihilates) an electron in a state i with spin a =f , J,, and energy £Vm) = e, =F Ez/2 
where is the orbital energy of i-th state. A <C 1 is the dimensionless interaction constant, and 5 is the average level 
spacing: 



ti+i 



(2.2) 



Let us stop for a moment to discuss the approximation made in Eq. (2.1). We included in Eq. (2.1) only the matrix 
elements of the interaction Hamiltonian responsible for the superconductivity. We omitted two kinds of diagonal 
terms. The term proportional to a iai a,i (T1 at %' CT2 represents the total charging energy responsible for the Coulomb 
blockaddlJ. It is not important for us because it does not lead to any anomalies at energies of the order of Zeeman 
splitting, and it can be accounted for by the corresponding shift of the applied bias. Other diagonal terms such as 
the one proportional to a\ ai a,i a2 aj CT2 a,j ai represents the spin exchange. It is not included because it leads only to the 
renormalization of the Lande factor g^. 

We also omitted off-diagonal terms, such as a^Oj <akiaii with i, j, k, I not equal pairwise, corresponding to the 
matrix elements: 

M% = J drdr'V(r ~ r')iP*(r)iP*(r')Mr)Mr')- 

The wave functions arc known to oscillate very fast, so the wavefunctions of different levels are very weakly correlated. 
We can restrict our consideration by short range interaction, V(r — r') = £-5[r — r'), where v§ is the bare DoS. 

One see that the integrand (|V't( r )l IV'j'WI ) m ^ ne diagonal matrix elements is always positive while the product 
(rp* (r)ip* (rtykirfipiir)) can be both positive and negative As a, result, the off-diagonal matrix elements turn out to 
be smaller than diagonal ones. Straightforward calculationl!3ll3 shows that they are smaller by the factor 1/g. 

In the paramagnetic state (Ez > \/2A), the structure of the ground state is similar to that without interaction, 
see Fig. ||. The orbitals with < —Ez/2 are doubly occupied while those with > Ez/2 are empty. The orbitals 
with |e,j| < Ez/2 are spin-polarized with up-spin. 
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FIG. 2. Structure of the ground state of the superconduc- 
tor above paramagnetic limit. Electron tunneling onto the 
orbital eo, creates the spin single states on this orbital. At 
some energy eo mixing of this singlet with the empty states 
becomes resonant, see the text. 



The Hamiltonian (2.1) does not affect the spin polarized states, but mixes the doubly-occupied and empty states. 
Since those states are separated from each other by a large gap E Zl this mixing can be treated perturbatively. Thus, 
the mixing does not change the ground state qualitatively. On the contrary the spectrum of the excitations, i.e., 
the tunneling DoS changes drastically due to the interaction. The essence of this effect is that spin-down electron 
tunneling into some orbital eo already occupied by spin down electron creates an electron pair which can mix with the 
empty orbitals and thus interact with superconducting fluctuations. This mixing turns out to be resonant at some 
energy E = E* and it leads to the sharp singularity in the spectrum of one-electron excitation. 

To evaluate the effect of superconducting fluctuations on the DoS of electrons in the paramagnetic state, we use 
the diagrammatic technique for the Green function (GF) at zero-temperatureO. The DoS can be expressed through 
the one particle GF, G i(T (cj), of an electron on the orbital j and with spin a = ±1 =f (J,): 

z/o-(u>) = -~sgn(Lj)lm.y^Gia(Lj), (2.3) 

TV * ^ 



where 



= Gl 1 - (2.4) 



G\ a is the GF for the non-interacting system 

G° m) = (u+ - ei ± E z /2)~\ (2.5) 

and Ej CT is the one particle self-energy. 

Leading contribution to the self-energy is shown in Fig. ||a. The solid and curly lines denote single-particle GF and 
the propagator of superconducting fluctuations, respectively. The latter can be obtained by summing the polarization 
loops in the Cooper channel shown in Fig. ||b. The single loop is given by 

1 , / _c^ 2 

where c<j+ = lu + iOsgn(uj) and uo c is the high-energy cut-off. Solving the Dyson equation (Fig. [|b), we obtain the 
propagator 



nM ^ ln Ur^i'' (2 - 6) 
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A(w) = 



XS 



28 



i - xsh(uj) ln (£1 



pa _,.,a \ ' 



(2.7) 



where A = cu c exp (— 1/A) is the BCS gap. 

The propagator (2.7) has the pole at w = ±fi, 



= ^|-A 2 . 

This pole can be interpreted as the bound state of two quasiparticles with energy f2. 

til + till 



(2.8) 




(6) 



'Tftftftr* = ■ + ^ n(u) \c5SCS^ 



FIG. 3. Diagrams for (a) self-energy £^1(0;) and (b) su- 
perconducting propagator A(a>). 

Analytic expression for the self-energy given by diagram Fig. ^a, has the form 

du>x 



2tt 



A( w + ^)G° T ( Wl ). 



(2.9) 



One can see that there are two contributions to the self-energy. One comes from the pole of A and the other is due 
to the branch-cut of this propagator. The pole contribution gives a singularity of the self-energy at certain to and ti 
while the contribution of the branch-cut is smooth. To find the singularity in the DoS, only the pole contribution to 
E may be retained: 



Ej|M 



<5A 2 



f 



VL uj + + e l - E z /2 + 0,sgn(ei - E z /2)' 



(2.10) 



At certain u> the pole of the self-energy coincide s wit h the pole of G°. This causes the singularity in the DoS. One 
can check that singularities of Eq. (2.1) and Eq. ( 2.10D coincide provided = 0/2 and 



E Z + n E Z + y/Ejr - A 2 _ 

= = tL 



Substituting Eq. (2.10) into Eq. (2.4) we obtain the GF for the down-spin electron at u> close to E*: 

r . s cj + + e i: - E z /2 - n 

l[UJ) ~ [lo+ -€i- E z /2)(w+ + e, - Ez/2 - n) - w$ ' 

where energy scale of the singularity is given by 



(2.11) 



(2.12) 
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w a = 



5A 2 



(2.13) 



Since Ez, A>i, one can neglect the fine structure of the DoS on the scale of 5 and substitute the summation over 
i by the integration over ef. 



E G u = ^ / 



<k. 



uj+ + e t - E z /2 - Q/2 
-e 2 + (c^+ - Ez/2 - n/2) - W$ 

to-E* 



-IVqTT 



- E*Y - W 2 



(2.14) 



Analogously, the GF for the up-spin electron can be obtained by changing the signs of Ez and f2, so that the 
singularity occurs at lo = — E* . 



Substituting Eq. (2.14) into Eq. ( |2.3| ), we obtain the final expression for the tunneling DoS in the ultra-small grains 

„2 \ 1/2 



. w ± E" \ I x- 

^T(;)M = t/ o J? o ( Wo J, F (x)=Re(^ 



a- 



1 



(2.15) 



wher e t^p is the bare Do S pe r one s pin, e nergy E* is defined by Eq. (1.7) and width of the singularity Wo is given by 
Eq. ( fTl3| ), Equations (§j|) and ( pl3| ) are the main result of this section. They predict the hard gap in the spin 
resolved density of states: v a {uj) vanishes at \u + aE*\ < Wo. Overall density of states + v-[ is suppressed by a 
factor of two near the singularity. 

In this calculation we neglected higher corrections to the self energy, e.g., those shown in Fig. ^a, b. In order to 
justify this approximation, wc ha ve to c ompa re contributions shown in Figs, ^a, b with the reducible diagram shown 
in Fig. [|c included in Eqs. ( 2.10 ) and ( 2.12 ). Singular contribution originates from the pole of A. It means that A 
carries frequency f2. The singularity in the DoS at ui — E* appears when the pole of the self-energy and the pole 
of G° coincide. This happens when GF G° for up-spin carries energy fl — to. In diagrams Fig. |^a, b, and c, the 
intermediate G° for down-spin should carry energy lo to give singularity to the DoS at E* . This condition can not be 
satisfied for the diagrams Fig. [|a, b. As a result, after the integration over the intermediate frequency, these higher 
order corrections turns out to be smaller than the reducible contribution (c) by small factor Wq/A ~ \J5jA~ -C 1. 




(<-) 




FIG. 4. Higher order correction to the self-energy (a) 
and (b), which were neglected in comparison with reducible 
diagram (c). 

III. DISORDERED INFINITE SYSTEMS 



In this section, we will obtain the quantitative results for the tunneling anomaly in the infinite systems. In subsection 



[II A wc will start from the perturbation theory and demonstrate that the lowest order perturbative results diverge 
algebraical ly at energies close to E* . In order to deal with this divergence, we develop non-perturbative approach in 
subsection [II B and obtain the analytic expressions for the shape of the singularities for all interesting cases. This 
machinery will be also used later in Sec. 
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A. Perturbative results 



The analysis of 0-D system presented in the previous subsection, is not directly applicable to superconducting wires 
and films because one can not approximate interaction Hamiltonian 



/oo 
dra\(r)a\(r)ai(r)af(r). 
-OO 



(3.1) 



by its diagonal matrix elements. Despite this complication, we will still be able to show that the singularity persists 
and remains at the same bias as that in 0-D. 

To describe this singularity we once again have to evaluate the effect of the superconducting fluctuations on the 
GF of electrons. First of all, we need to evaluate the propagator of superconducting fluctuations A, see Fig. ||. In 
contrast with 0-D case, the superconducting fluctuations in the bulk system are inhomogeneous: the propagator for 
the superconducting fluctuations depends on the wave vector Q. (We will omit the vector notation in the momenta, 
e.g., Q = Q.) Solving the Dyson equation, shown in Fig. Ha, we obtain the propagator 



A(w,Q) 



z/ rn ( 



£j-(M+i-DQ 2 ) 5 



(3.2) 



where D is the diffusion coefficient. At Q = 0, propagator (3^2) resembles zero-dimensional expression (2/7). We see 
that the propagator has the singularity at to close to ft provided DQ 2 -C O. As we will see shortly , it r esults in the 
singularity in the DoS developing at exactly the same energy as in 0-D case, to — E* , given by Eq. ( 2.11 ). 



A (w, Q) 



Mr 
□ 



P^^a5OTTi (a) 




('') 



: %G m (u, P ) : 



1 

2lTV T 



uj - e p ± E z /2 + sgn(ui) ^ ' 

FIG. 5. Diagrams for the (a) propagator of superconduct- 
ing fluctuations A (ui, Q) and (b) the vertex function j+(u>, Q). 

Next step is to consider vertex function in the particle-particle channel. The ladder approximation which gives the 
main contribution at epr 3> I ( t is the elastic mean free time) is shown in Fig. ||b. Analytically, they are given by 



Here I±(lo, Q) stands for 

/±(w,Q) = 



f 



I 



7±(w,Q) = t + I±(ui,Q)j±(ui,Q), 

(*)Gt(i)(w + wi, e p )G im (u)i, e_ p _ Q ) 
(dp)- 



(3.3) 



2m> T J y rl (u + uji - e p ± + -^sgn(uj + uj 1 ))(lo 1 - e_ p+Q + ^sgnluj^) 

= {l + T [i(\uj\±sgn(Lj)Ez)-DQ 2 ]}6[-(u J + u Jl )u J }, (3.4) 

where (dp) — , and the GFs are averaged over disordcrSR Here we used the conditions of the diffusion approxi- 



mation, lot <C I and Ql <C 1. We substitute Eq. (3^) into Eq. fl3.3| ) and we obtain 

6[- (w + wi)u] 



j±(oj, Q) = t6 [(uj + uji)oj] + 



-i[\w\ ±sgn(uj)E z ] + DQ 2 ' 



(3.5) 
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FIG. 6. Diagrams for the lowest order corrections to the 
one-particle Green function. 

Now we can evaluate the first order correction to the one-particle GF, Fig. ||: 

SGf\u,p) = -^ J(dQ) J ^Gl(cj,p)A(cj + cj 1 ,Q)G r (u; 1 ,-p + Q) 1 2 _(u 1 -uj 1 ,Q) 

The DoS is determined by the GF integrated over p: 

sgn(uj) 
"TU)( W ) = ; 



Im / (dp)G ul) (w,p) 



(3.6) 



(3.7) 



We substitute Eq. ( |3.5| ) into Eq. (3/7) and perform the integration over p and u>\. For Ql -C 1 and \u> — wi|t *C 1 
integration over the momentum p results in 

(dp)G^(u,p)Gi(oJi > — p + Q) = i2TTVQ6(—LOLUi)T 2 sgn(uji). 

Performing integration over ui\ we take into account only the pole contributions in the propagator ([T^) for DQ 2 <C fl: 

1 1 



A 2 

A(c,Q)« — 



v n \n + uj-iDQ 2 Vl-uj-iDQ 2 



(3.8) 



since the integrals along the branch cuts give only the corrections which are smooth functions of u). The main 
cont ribu tion to the frequency integral results from the region where the real part of the pole of the propagator A in 
Eq. (3.7), Recji = — lo ± f2 is close to that of the vertex 7_, Re lo\ = lu + sgn{u>i — lo)Ez, and the imaginary parts of 
those poles have different signs. The latter requires u> R ecjj < 0, us 2 > [Keuji] 2 . One can easily check that all these 
conditions can be met only if u> is close to E* from Eq. (1.7). 

Evaluating the integral over uj\ in Eq. ([3.7]) we obtain the first order correction to the DoS 



M;{)H _ a 2 

vq 2i/qQ 
where Cjq)(w, Q) is the Cooperon given by 

C(oj,Q) = 



Re J{dQ)C 2 {Lo±E*,Q), 



DO 2 



(3.9) 



(3.10) 



[Calculation of 8v^ requires an obvious modification of Eq. ( |3.7[ ).] 
For one dimensional case (wire), this correction acquires the form 



Sis 



(i) 
T(-l) 



A 



3/2 



(3.11) 



It is possible to neglect higher order corrections to DoS only provided |w±i?*| is large. For uj — > E*, correction (3.11) 
diverges. Therefore we need to sum up all the orders of perturbation theory to describe the DoS in the vicinity of E* . 
Such calculation is carried out in the next subsection. 
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For two dimensions (films), the first order correction to the DoS vanishes for |u;| 7^ E*. However, this is nothing 
but an artifact of the first order approximation and the second order correction is already finite. Diagrams for this 
correction are shown in Fig. [?]. The result can be written as: 



8vW{w) ul) A 2 \ 2 V 



= -2 . 

For two dimensional case Eq. ( p. 12 ) gives 



— Im / (dQi)(dQ 2 )C' m) (uj,Qi)C m) {cj,Q 2 ). 



ou Ul) 



hi 



n 



\lo±E*\ 



(3.12) 



(3.13) 



where g = AttDuo 3> 1 is the dimensionless conductance of the film in the normal state. Deriving Eq. ( 3. IS ) , we cut 
off the logarithmic divergence a t la rge momenta Q2 by the condition DQ\ < A, since it determines the applicability 
of a single pole approximation (3 






p + Q,u> + wi,t 




— i sgn(u)) E z + DQ' 2 



FIG. 7. Diagrams for the second order corrections to the DoS. Diagrams irreducible with respect to the curly lines, similar 
to Fig. ^|a,b, are negligible for the reason discussed in the end of Sec. [n|. 

As well as II? case, the perturbation theory fails t o de scribe the DoS in the vicinity of E* in two dimensions. It 
is noteworthy that the singularity described by Eq. (3.13) is much more pronounced than that of the normal metal 
(T > T c ) which arises due to the superconducting fluctuations and it is of the order of g _1 ln[ln(a; ± Ez)}, see Rcf. H. 
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The enhancement of th is si ngularity results directly from the isolated pole in the propagator of the superconducting 
fluctuations, see Eqs. (3.2) and (ETa) rather than in the branch cuts of Ref. ||. 



B. Non-perturbative results 

1. Derivation of self- consistency equations 

We start summation of the perturbation theory terms from the series of diagrams for DoS presented on Fig. ||. 
Such diagrams dominate in each order of the perturbation theory in comparison with those of the same order but 
including irreducible with respect to curly lines elements, similar to Fig. ^a,b. The later statement can be justified 
by arguments similar to the analysis of diagrams Fig. || presented in the end of Sec. O. 




FIG. 8. Structure of the fc-th order perturbation theory. 
Vertical wiggly line on the upper Green function corresponds 
to cutting this GF into two lines and fixing its frequency u>. 



The diagram of the k-th order of the perturbation theory, Fig. [|, contains k curly lines that stand for the fluctuation 
propagator K(ujj,Qj) in Eq (3.2) (1 < j < k). It also contains (2k + 1) GF. Before averaging over disorder, each 
GF G a (ijj;p,p') depends on two momenta (initial p and final p') and on the direction of the spins a = ± =t,|- 
Contribution of this diagram to the DoS of j electrons 8vi(lo) at positive energy u> > is 

= -^Im / J {dpi)Wl) . . . {dpk+l) J (d Ql )(dQ 2 ) ■ ■ ■ (d Qk ) 

k 

x Gi(uj]p k+ i,pi) Y[ A(uj,Qj)Gi(uj;pj,Pj)Gj(-uj + ujj; -p'j + Qj, -pj+i + Qj). (3.14) 



Disorder averaging of the product of the Green functions in Eq. ( [3.14 ) and of the superconducting propagator A 
can be carried out independently, since their correlation gives rise to non-singular corrections containing additional 
smallness 1/(£ft)- 

For the sake of convenience we introduce G R ^ A \uj]p,p'): retar ded ( advanced) GF at Ez = 0. In the absence of 
spin-orbit scattering the Green functions G R ( A \ui]p,p') from Eq. ( 3.14 ) can be presented through G R ^ as 



G a (uj;p,p') = 



(G% 
\G A ( 



^ + aE z /2;p,p') uj > 
^ + aE z /2;p,p') uj < 



(3.15) 



31 



After substituting Eqs. j3.15| ) and (|J) into Eq. (fy]i|), one can integrate over all intermediate frequencies uj 
because the fl uctuation propagator A (a;,, Qj) has simple poles at ujj = ztfl^fiDQ^. According to Eq. (3.15), G[(ui;p,p') 
in Eq. ( 3.14 ) is retarded (ui > 0). In order to get the advanced GF from Gf (uij — ui; —p'j +Qj 7 —pj+i + Qj) we choose 
the pole in A(wj-; Qj) with positive real part, uj 3 ■ — fl — iDQ?. (One can check that another pole leads to the product 
of the retarded functions, which vanish after the disorder averaging.) Using Eqs. (3.8) and (2.8) and introducing 
short-hand notation 
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-u + 2E* 



iDQ 2 



(3.16) 



present 5v^(uj) at uj > as 



(. 



1 



7r \ i/ r2 



~ n 

Im / (^ +1 )G fl (w; Ffc+l!Pl ) J](*i)(^)Wi) 



G^^^G 4 ^; + Qi, + Q 3 ). 



Averaging Eq. (3.17) over disorder and using the identity 



(dp 1 )G R (u;p 1 ,p' 1 )G R (u;p k+1 ,p 1 ) = -—G R (uj;p k+1 ,p[), 



(3.17) 



(3.18) 



we obtain 



1 / A 2 



7T \^0^ 



1 



d 



2k J duo 



■Im(/nWi)fe)(^: 

i=i 

x G R {uj;p h p' j )G A {ujQ j ; -p' } + Q 3 -, -p 3+ i + Q 3 -)). 



(3.19) 



The role of the factor 1/fc in Eq. ( |3.19] ) is to cancel the multiple counting: there are k retarded GF in Eq. (3.19), 
and application of the operator d u to any one of them leads to Eq. (3.17). In addition, Eq. ( 3.19D includes terms like 
dG A /duj. One can check that they give exactly the same contribution as terms which contain dG R /duj. Additional 
factor 1/2 takes these contributions into account. 

Using a conventional trick 1/k — rj k dr], (see Ref. [7]), one can present (Sv^ (ui)) in a form 



-LA-I m jL / dr] I (dQ)(dp)(G H (u!- 1 p)) (G a (ujq; —p + Q)) 



2ttt vq^I duj J 

x 7i fc) (^Q,Q), 



where short hand notation 



jq = oj - uj q = 2lu - 2E* + iDQ 2 



(3.20) 



(3.21) 



is used and ujq is defined in Eq. (3.16). Similar equation holds for spin up and negative to. In this case one should 
use ujq = -2oj - 2E* + iDQ 2 . 

The vertex can be written using Fig. |^. The rules of reading diagrams on Figs. ^ and ^ are slightly different 
from the conventional rules we used before: (i) to each curly line corresponds factor r]A 2 / (i/qCI); (ii) no summation 
over the f requencies is implied: each retarded GF bears frequency u, each advanced GF bears frequency ujq given by 
Eq. ( 3.16| ); (iii) each interaction with curly lines changes retarded Green function to advanced and back. The resulting 
expression reads: 



jW(u> q ,Q)=t(t 1 



A 2 \ 



k-l 



k-1 



HidQMdp^dp'j) 



G R (uj;p j ,p' j )G A (uj Qj ; -pj + Qj,-p j+ i + Qj) 
x G r (uj;p' 1 , P )G a (ujq; -p' k + Q, -p + Q)). 



(3.22) 



Note that the averaged product of GFs in the integ rand of Eq. (3.22) does not depend on p. We can thus perform in 
Eq. ( 3.20 ) the integration over p. As a result, Eq. ( |3.22 ) takes a form 



— = -— 7^ lm S- [ df] I (dQ)j(u)q,Q), 

l/Q VQil OUJ Jo 



(3.23) 



where 



K^Q,Q) = ^7 (fc) (^Q,Q)- 



(3.24) 



fc=i 
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FIG. 9. Structure of the 7^. Three points defining vertex 
are denoted by thick dots. 

Our goal is to evaluate 7(£>q, Q) self-consistently. The difficulty is that diagrammatic series for 7 contains other 
elements except 7 itself, see Fig. [LQ, 




FIG. 10. Diagrammatic equation for the vertex "fu>((,Q)- 
These elements are known as Hikami boxesEl. The simplest Hikami box B^ 2 \ which appears already in second 



order of the perturbation theory for 5v, is the integral (dp) of a sum of three diagrams shown on Fig. 1 1 

1 



B< 2 >(wi,w 2 ,Qi,Q 2 ) 



(2w ) 3 



-i(oJt+u 2 )+D(Qi + Qi)}. 



(3.25) 



(2 I ti/t) 4 B( 2 1(u 1 ,u 2 ,Q 1 ,Q ;! ) 




V X 

X X, 



Q K {u,p) \ / G A ^-^ 2 ,-p+Q 2 ) 



XX + XX 

xV + xV 
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FIG. 11. Second order Hikami box B (2) . 
The third order Hikami box, B^ 3 \ is given by the momentum integral (dp) of a sum of sixteen diagrams of Fig. nl 



B (3) (wi,w 2 ,w 3 ;Qi,Q2,<23 



(2m/„) ; 



: H(cji + u 2 + lo 3 ) + D(Qf + Q 2 2 + Ql)]. 



(3.26) 




V V 
A A 



X 



\ 
/ 



. 5 rotations of the 
dashed line 



. 2 rotations of the 
dashed lines 



// X 

x : -/ 



5 rotations of the 
dashed lines 



FIG. 12. Third order Hikami box B (3) . Momenta and 
frequencies in the Green functions are arranged similar to 
Fig. 0. 



The equation for the vertex j(ujq,Q) = 7 W (C = &Q,Q) that determines correction to the DoS [Eq. ( [3.23 )] 



7<-(C. Q) = i + - + K - D Q 7«(C. Q) 



V ) I Wi)(2^o) 3 S (2) (^Q^Q 1 ,Q,Qi)7 2 (l)7-(C,Q) 



'/"( — ] / / (dQ 1 )(dQ 2 )(2Trv ) 5 B( 3 \u Q ,u Ql ,u Q2 ,Q,Q 1 ,Q 2 ) 



\vaSiJ 

x 7 2 (l) 7 2 (2) 7w (C,Q) 



(3.27) 



where "f(j) = "fui^Q^Qj) with CjQ j determined by Eq. ( 3.21 ). We introduced extra variable £ in order to separate 
the external energy u> and the integration variable in Eq. (3.27) and further. 



Hikami box of the k-th order B^ is given by a sum of diagrams (of the type of Figs. 11 and |l2|) which contain 2k 
vertices. Strictly speaking, B^ depends on 2k — 1 sets of momentum and energy transfer (uj,Q). However, in order 
to evaluate the DoS correction given by diagrams Fig. 0, we can restrict ourselves by Hikami boxes that depend only 
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on k sets (u>j, Qj), 1 < j < k: (u>j, Qj) and (— ujj, —Qj) characterize neighboring vertices. Equations (3.25) and ( 3.2(i| ) 
allow us to conjecture that B^ at arbitrary k has a form 



c k 



(2™ ) ; 



(3.28) 



where Ck are numerical coefficients. 

We are not going to determine coefficients Cf. directly. An important feature of Hikami boxes is that they are local 
obje cts, d etermined by distances smaller or of the order of the mean free path I. Therefore, the coefficients Ck in 
Eq. ( 3.28) ) do not depend on the dimensionality. Therefore, we can compare the exact solution of 0-D problem, Sec. 0, 
with the su m of p erturbation theory series involving coefficients Ck and thus find those coefficients. 

Equation $.2% can be rewritten in terms of Ck as 



00 r A 2 f 

He + dq 2 ) 7u (c, Q) = i + HC + DQ 2 hUc, Q) J2 c * [v^n J WOt(i) 

00 r a 2 f ~\ k ~~ 

Introducing the function, f(x) 

oo 
k=l 

we obtain a simple equation for 7^(C,Q) 

(-»C + DQ 2 ) lu] (C, QOTo) + 7«(C, Q)ft/'(A>) 
where /3q and /3i are connected with 7 by 

A 
A 



(3.29) 

(dQ 2 )(-iu Q2 +DQl) lw {u Q2 ,Q 2 f. 

(3.30) 

1, (3.31) 



/?i(w) = i] 



{dQ^iiCjQ^^-iuQ + DQ 2 



One can use Eq. ( |3.3l| ) to present 7^(C, Q) in the form 

7-(C.Q) = 



Z(w) 



-< + £>Q 2 + 2m(w) ' 



where parameters Z and m are determined as 

1 



/GW) 



, 2m(u)=A(w) 



/(AH) ' 



(3.32) 



(3.33) 



(3.34) 



It was already mentioned that the Hikami boxes B^ and thus coefficients C n as well as the function fix) do 
not depend on the number of dimensions d. In contrast, w-dependence of parameters (3q, Pi and m are different at 
different dimensions. Let us first consider 0D case in order to determine f{x) explicitly. At d — one has to abolish 
integration over Q in Eq. ( 3.32 ) and substitute inverse mean level spacing <5 _1 for vq and 2(u> — E*) for Qq. The 
vertex 7 can be determined straightforwardly 



7.(0 



(3.35) 



wher e 0D GF Gi a and G\ a are determined by Eq. (|2.12|) [Eq. (2.13 ) for Wq should be multiplied by y/rj] and by 
Eq. (2.5) respectively. Substitution of Eqs. ( 2.12 ) and (|2.5D into Eq. ( 3.35 ) gives after the integration over e. 



7.(0 = g 



uj-E* 
W 



-iC + i( w -E*)[l-l/F (^)] : 



(3.36) 
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wher e the function Fp(x) is given by Eq. (2.15), E* is determined by Eq. (1.7), and Wq — r]5A 2 /ft. By comparing 
Eq. ( 3.36 ) with Eq. ( 3.35 ) we immediately obtain Z(u>) and m(w) for zero dimensions: 



F { 



j-E* 
W„ 



Z q (cj) = 



mo(w) = i(uj — E*) 



On the other hand, from zero dimensional version of Eq. (3.32) at £ — 2lu — IE* we can determine {3q and f3\: 



(3.37) 
(3.38) 

(3.39) 



N ow we are in state to determine the functi on f (x) from Eq. ( |3.30| ). We express Fq through /3q, substitute it into 
Eq. (3.37) for Zq and use the connection Eq. (3.34) between Zq(lu) and f(0o). As a result we have 



2x F (2x) 



1 

2x 



yi -Ax 

2x 



(3.40) 



This function al de pe ndenc e w hich remai ns the same at all dimensions, can be used to evaluate 6v±(u>) for d = 1, 2. 

Equations ( 3.23| ), ( |3 . 32 ) - (3.34 ) and ( 3.4C ) constitute complete set allowing to find the DoS in any dimensions. 
One has to substitute Eq. ( 3.33| ) into Eq. (3.32), and find functions Z(u>) and m(u>) self-consistently with the help of 
Eqs. ( 3.34 ) and ( 3.4C ). The result should be substituted in Eq. ( 3.23 ) which gives the final non-perturbative answer 
for the DoS. In the following subsection, this program will be carried out for ID (wires) and 2D (films) systems. 



2. Solution of self- consistency equations 



We substitute Eq. ( 3.33 ) into Eqs. ( 3.32 ) and ( 3.23 ) and perform integration over the wavevectors Q. Equation 
( |3.23 ) acquires the form 



Here, we used Eq. (3.34) and introduced dimensionless frequency and mass 

u> + <jE* m 



W d W d 

The relevant energy scales, which, as we will see below, are the widths of the tunneling anomaly, are given by 

v2 X 2/3 A 2 



Wi = 3 



(- 



A 2 



\l6vinVD 



w 2 



(3.41) 



(3.42) 



(3.43) 



where is given by Eq. (2.8), D is the diffusion coefficient, v\ — (mp f S) / {2ir 2 ) is the one dimensiopal density of 
states per unit spin (S is the cross-section of the wire), and g — Att^D is the dimensionless conductance^. The latter 
is related to the normal state resis tance of the film as g = 25.4fc£l/i? a . 
Dimensionless functions in Eq. fl3.41 ), Md(x), are defined as 



Mi (a;) = 



, M 2 (x)=ln(^ 



(3.44) 



To find M2 we cut off the logarithmic divergence at la rge momenta Q by the condition DQ 2 < A, which determines 
the applicability of a single pole approximation (3.8), and neglected the factor A/51 ~ 1 in the argument of the 
logarithm. 



Using the same notation, we obtain from Eqs. ( 3.32 ) and ( 3.34 ) 

/3o(^){/[/?o(^)]} 2 = - V M' d [-iu; + m(uj,r,)} 
m {f\po(w, v)}} 3 = vf'Woi^ V)] Wd + m(uj, V )} + m{uj, r})M' d [-iu + m(co, V )}} 



(3.45a) 
(3.45b) 
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Equation ( 3.40| ) allows to formally solve Eq. (3.45a): 

r)M' d (—iuj + to) 



Po = ■ 

[1 - rjM d {-iu + m)Y 

f(fa) = l-riM' i {-iw + m), 

f'Wo) 1 



[/(/?o)] 3 l + vM' d (-iu + m) 



We can now substitute Eq. ( 3.46c ) into Eq. ( 3.45b ) and obtain after simple algebra 

m(u>, rf) = r\M d [—iuj + m(w, rj)} . 



(3.46a) 
(3.46b) 
(3.46c) 

(3.47) 



Further calculation are substantially simplified due to the fact that the integrand in Eq. (3.41) can b e pres ented 
as a total derivative with respect to r\. In order to demonstrate this, we differentiate both sides of Eq. ( 3.47 ) with 
respect to r\ : 



dm dm 
— — = M d (-iu; + to) + r]M d (-iu + m) — — . 
or) or] 



(3.48) 



Find ing d m/drj from Eq. ( 3.48j) , we notice with the help of Eq. ( 3.46b| ) that it coincides with the integrand in 
Eq. fl3.4l| ). Integration in Eq. (^.41|) can be immediately performed and we obtain 



5v a 



~2— — Im m(u!, rj 
ouj 



(3.49) 



Finally we put rj = 1 in Eq. (3.47), differentiate both sides of this equation with respect to to, and substitute the 
result for dm/duj into Eq. ( 3.49| ). After restoration of original units for u according to Eq. ( |3.42 ), we obtain for the 
density of states v a (uS) = vq + Si/ a (u>) the following result 



v„{w) 



F d 



u) + oE* 

Wd 



(3.50) 



where a = ±1 corresponds to th e sp in-up and spin-down densities of states respectively, and the widths of the 
singularity W d are defined in Eq. ( 3.43| ). Dimensionless function F d {x) is given by 



F d {x) = Re 



1 + z(x) 



1 — z(x) 

where function z(x) is implicitly defined as the solution of equations 

z{x) = M' d [—ix + y(x)] , 
y(x) = M d [-ix + y(x)] , 



(3.51a) 



(3.51b) 



with functions M d (x) being defined in Eq. (3.44). If Eq. (3.51b) has several solutions, one has to choose the one 
reproducing the pe rturba tion theory at i > 1 and remaining on the same list of the Riemann surface at small x. 

In ID case, Eq. (|3.51h| ) can be rewritten as a cubic equation and solved using the Cardano formula which yields 
universal (independent on and D) expression for the shape of the singularity: 



Fx {x) = l- ^Rc [l - ix ( Vl + y(x) + tyl - y(xj) 
y{x) = vl + ix 3 . 



Functions ^fz in Eq. (3.52) are defined to map compl ex p lane — ir < a rq(z) < ir to — ? < arg(z) < ~. 
For the two-dimensional case, we obtain from Eqs. ( 3.51 ) and ( 3.44 ): 



F 2 (x) = Re 



1 — z(x) 
1 + z(x)'' 



(3.52) 



(3.53) 



where z(x) is the solution of the transcendental equation 
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z{x) = 



1 

-ix + ln[4pz(a;)] 



(3.54) 



Shape of the singularity in this case depends on the cond uctan ce g and, thus, is not universal. However , thi s 
depe nden ce is rather weak. For \u> ± E*\ Wd, Eqs. ( 3.5 2| ) and (3.53) match the perturbative results, Eqs. (3.11) 
and ( 3.13 ). Found energy dependence of the density of states is shown in Fig. [l3| 








I 








-0.2 


















v 100 /// 






-0.4 


















\ 50 / 












\ 20 / 


(b) 




-0.6 






VVj) = 5 








-10 


-5 


5 


in 


ii 



FIG. 13. Singularity in DoS for spin-down polarized elec- 
trons for (a) ID, and (b) 2D syst ems. Widths of the singu- 
larity W 12 ar e given by E q. (3^43 ), and the shape is defined 
by Eqs. (po|), (isjf) and (3.53[). 



IV. QUALITATIVE DISCUSSION 



In this section we will present a simple qualitative interpretation of the main results obtained in the previous 
sections. We believe that this simplified way of thinking provides instructive physical intuition even though it fails to 
give completely quantitative description. 

It has been already emphasized in the beginning of section Sec. |0| that the ground state of the system above the 
paramagnetic limit has the structure similar to that of a non-interacting system. All the mixing of the noninteracting 
states caused by the interaction part of the Hamiltonian ( |2.1[ ) is perturbative. We can neglect it completely in a 
rough approximation and consider the electrons occupying orbitals with orbital energies < —Ez/2, see Fig. ^ to 
be "frozen". As soon as the spin down electron tunnels onto an orbital < eo < Ez/2, see Fig. the electron pan- 
on this orbital is created. Due to the interaction, this pair can mix with all the empty states, ti > Ez/2. It is this 
mixing which gives rise to the singularity in the DoS. On the other hand, within the same approximation, all the 
electrons on the orbitals e; < —Eal2 can still be considered as "frozen" . ( This approximation is similar in spirit to 
the well-known Cooper proceduretEl. ) 

Thus, we arrive to the following recipe for the evaluation of the energy of one electron excitation. First we have to 
find the eigenenergies E-L-. (eo) of the two-electron problem within the Hilbert space consisting of orbital eg and of all 
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the orbitals e, > Ez/2, see Fig. [l4|. (This energy spectrum naturally depends on e as parameter.) Then, the energies 
E J , (eo) of the one particle excitation, corresponding to the introduction of electron onto the orbital eo are 



£|(eo)=^ 2) (eo) 



Ez 
2 



(4.1) 



since the total energy of the electron which occupied this orbital before the tunneling event was eo — Ez/2, while the 
state of the rest of the electrons did not change. Accordingly, the density of states for spin-down electrons is given by 



(4.2) 







FIG. 14. Hilbert space for the solution of two - electron 
problem. All the orbitals (dotted lines) with < —Ez/2 are 
excluded since they are occupied by "frozen" electron pairs 
and the orbitals —Ez/2 < e; < Ez/2,et / eo are excluded 
becaus e sin gle occupied orbitals are not affected by the inter- 
action (2.1). 



Now, we have to find the spectrum of two electron eigenenergies E J ^(eo). Since the interaction in the Hamiltonian 
involves only the spin-singlet orbitals, the wave function of the electron pair ip can be labeled by one orbital index 
and it is governed by the Schrodinger equation: 



E^ipi = 2e i '4> i - X6^2 i>j 



(4.3) 



The eigenenergies E J ,^ can be determined from the following equation 



2e n - E 



(2) 



E 



2e- - FJ 

e z >E z /2 -°(2) 



(4.4) 



For low-lying eigenstates E 3 ,^. < Ez, one can substitute the summation in Eq. (4^4) by the integration. Given the 
high-energy cut-off uj c , it yields 
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where Af, = oj c exp (— 2/A) is the binding energy of the Cooper pair. 



As we will see in a moment, 



2e 



E 3 

^(2) 



it can be expanded in the Taylor series, ln(l 



3> 8, and, therefore, the logarithm in Eq. (4.5) should be also small and 
x)k 
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\x\ <C 1. Equation ( 4^5 ) is immediately simplified to 

E z -A b 



^(2) 



2e 



^(2) 



A, 



and we obtain the solution for two relevant eigenenergies: 



Ef 2) (e a ) 




(4.6) 



All the other two-electron states have energies larger than Ez and they are not important for us. In Eq. (4.6) energy 

n b = E z - A b (4.7) 

has the meaning of the energy of the bound state of the Cooper pair measured from the Fermi level. It plays the 
role of e nergy fr om Eq. (2.8) in the rigorous solution. We will return to the discussion of the discrepancy between 
Eq. (fh7|) and Eq. (pit) later. 

Substituting Eq. ( |4.6|) into Eq. (4T), we obtain the energy of one particle excitations 



Ef(e ) = E b *± 




25A b , 



where the position of the singularity 



Ez + ttb 



(4.8) 



(4.9) 



is similar to the energy E* in exact Eq. ( [2.11 ) up to the substitution Q ^ fl b . 

According to Eqs. (4.2) and Eq. (fO|), the density of states van ishes in the energy strip \E — E£\ < (25Af ) ) 1 ^ 2 - 
hard gap in the DoS is formed, compare with the exact result, Eq. ( 2.15| ). The origin of this tunneling anomaly is the 
avoided crossing of the two-electron state formed by the tunneling electron (energy 2eo) with the bound state of the 
Cooper pair (energy fl b ). 

It is also noteworthy that if a spin-up electron tunnels into the grain, it never finds the pair for itself, and, 
therefore, no tunneling anomaly happens in this case. It means that the overall DoS does not vanish but rather shows 
the suppression by a factor of two. However, for the spin-polarized electrons tunneling into the grain, we predict the 
complete suppression of the tunneling DoS. 

The same arguments allow to justify the similar singularity, when spin-up electron with energy —Ez/2 < E < 
tunnels out from the system, while the spin down electrons tunneling from the system are not affected. 

The qualitative consideration above grasps the correct physics, however it fails to describe the effect quantitatively, 
it predicts correctly neither the position nor the width of the gap. This is similar to the discrepancy between the 
binding energy Af, in the original Cooper procedure and the correct BCS gap A: all the electrons below the Fermi 
energy were frozen. To remedy this drawback, one has to employ a parametrically exact procedure described in Sees, 
and 







III 



Let us now discuss the results obtained for the disordered bulk systems, Sec. [II. They can be briefly summarized 
as follows: (i) Singularity in the DoS persists; (ii) its posit ion does not change; (iii) energy scale of the singularity 
depends on both disorder and dimensionality, see Eq. (3.43). 

In order to understand the physics behind the singularities in the bulk systems, let us recall the meaning of zero- 
dimensional approximation. Strictly speaking, it implies that during the time tE — ^/E, (where E is the energy 
scale relevant for the problem), diffusively moving electron can visit all the system. The characteristic time of this 
diffusion is ~ L 2 /D, (L is the characteristic size of the system) which means that zero-dimensional approximation is 
applicable provided that the energy scale E does no t exc eed the Thouless energy E c — %D/L 2 . In our problem the 
relevant energy scale is the gap width Wo from Eq. (2.13), and condition 
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1 > W° oc L 2- d /2 

E c . 



(4.10) 



is definitely violated for the infinite systems L — > oo, (here d = 1, 2 is the dimensionality of the system). 



It is clear, that as soon as the condition (4.10) breaks down the geometrical size of the system L as well as its 
mean level spacing 5 = l/{vQL d ) is not relevant since the electron can not diffuse during finite time tjg = h/E over 
the distance larger than Le = \JDtE- On the other hand, it effectively visits all the space on the scale smaller 
than Le- Therefore, the following approximation holds: in order to understand the properties of the diffusive system 
associated with the energy scale E, we can separate the system into smaller patches of the size Le = \/TiD /E, and, 
then, apply zero-dimensional description to each patch independently, (assuming that different patches do not "talk" 
to each other). 

Let us no w apply this strategy to the problem in hand. First, we notice that the position of the singularity in 0D 
grain ( 2.11 ) does not depend on the size of the grain and therefore the singularity in each patch should be at the 
same ener gy E* as in zero-dimensional systems. Second, level spacing S entering into the width of the singularity 
Eq. ( 2.13 ) does depend on the size of the patch 



6 = 



voL d E 



(4.11) 



In this formula, scale Le is itself determined b y the width of t he sin gularity, E ~ Wd, so that the scale Wd has to be 
determined self-consistently. Substituting Eq. ( 4.11 ) into Eq. ( 2.13|) , we find 



( A ' 2 



K^oL^H 



1/2 



D 

w d 



1/2 



(4.12) 



Solving Eq. ( pLl^ ), we obtain 



A d/2 



flD d / 2 



(4.13) 



which agrees with the rigorous results ( 3.43 ) for one- (d = 1) and two- 
quantitative description requires machinery like one used in Sec. III. 



(d = 2) dimensional systems. However, the 



V. RELEVANT PERTURBATIONS 



A. Spin-orbit scattering 

In our previous consideration we assumed that electronic spin is a good quantum number, i.e., impurity scattering 
of electrons does not cause spin-flips. There are two sources of spin relaxation of conduction electrons: localized 
spins (paramagnetic impurities) and spin-orbit (SO) scattering of electrons by non-magnetic disorder. The latter is 
characterized by the scattering amplitude 

iv ao (]pfXPi]-tr)/p%, (5.1) 

where and p i are final and initial momenta of an electron, and cr is the spin operator cr = (cTxi&yT&z) whose 
components are Pauli matrices. It acts on the spinor wave function of the electron. 

Let us discuss the effect of SO scattering first starting with the qualitative physical picture. In the absence of both 
SO interactions and magnetic field two spin states which belong to a given orbital have the same energy. Magnetic 
fields splits this degeneracy. It was important for us that the splitting Ez is exactly the same for all of the orbital 
states. Now let us turn on SO interaction r _.Without external magnetic field the states are still double degenerate 
due to the T-invariance (Kramers doublets^). Magnetic field is well known to split the Kramers doublets similar 
to how it splits pure spin states in the absence of SO interactions. The main difference is that this splitting is not 
exactly uniform any more (see e.g. Ref. |l^). It is this dispersion of splittings that smears the DoS singularity. Zeeman 
splitting dominates the magnetic field effect on superconductivity only provided SO interaction is weak. However the 
DoS singularity turns out to be sensitive even to a weak SO scattering, since the characteristic SO energy(dispersion is 
splitting of Kramers doublets) should be compared with the width of the singularity Wd rather than wit h the s plitting 



Ez itself. The DoS for finite SO scattering can be evaluated in a way similar to our calculation in Sec. Ill B 
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Cooperon (or vertex) is formed by two electron Green functions. In the absence of external magnetic field it is 
convenient to classify Cooper poles by the total spin of two electrons S+ — {er\ + o"2)/2. Spin-orbit scattering does 
not affect the spin singlet part of the Cooperon {S 2 + = 0) because the spin-orbit scattering preserves T- invariance. 
However, this scattering leads to total spin relaxatioa, i.e., triplet (S 2 + = 2) component of Cooperon decays (pole in 
w-plane is shifted from the real axis even for Q = 0)E3. 

External magnetic field is coupled with the difference S- = (cri — er 2 )/2 of two electron spins, and as a result 
we classified Cooperon by the eigenvalue of the operator S- ■ Ez- These eigenvalues for S 2 _ — 2 are Ez, 0, — Ez 
corresponding to = 1, 0, —1 and for S 2 _ = 0. Neither of those two classifications is exact when both magnetic 
field and SO scattering take place simultaneously: operators S+ and S*L do not commute. On the other hand, as it 
was already mentioned, we should assume that SO effect is weak. This allows us to eval uate Cooperon perturbatively. 

Summing usual ladder diagrams and taking additional scattering amplitude Eq. ( |5.l| ) into account, we end up with 
an equation for 4x4 matrix of the Cooperon which we already discussed qualitatively: 



( 



DQ 2 )l + iE z -S- 



2S Z + 



6-'-, 

T 



(5.2) 



where t so = 1 / (2-kvqv 2 ) is the time of the spin relaxation by SO scattering, matrix I is the direct product / = <7g (g><7g, 
where unit matrix <7q is acting in the spin 2x2 space, and <7q is a unit matrix in 2 x 2 space of the electron lines. 
Operators S± are defined as 2S+ = (cr s + tictq) ® erg and 2S- = (cr s — tictq) ® <rf , where n = (l x , l y , l z ) is the unit 
vector, and <x s = {a s x , tr*, cr|) are the Pauli matrices in the spin space. 



Instead of diagonalizing Eq. (5.2), we can just evaluate correction to C 



HI)' 



c 



(5.3) 



in the first order of the perturbation theory in (Ezt so ) 1 < 1. This correction turns out to be 2/(3r so ). As a result, 
SO scattering simply shifts the Cooper poles C^ny. 
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-iu T iE z + DQ 2 + 



(5.4) 



From Eq. ( ^.4| ) one can guess that the all the interesting results can be obtained from the results of Sees. || and III 
by substituting 



iT, 



(5.5) 



where T is the spin orbit rate 



3t,, 



(5.6) 



so that the final result for the density of states is 



y T(l) 



RcF d 



uj±E* 



iT 



W d 



(5.7) 



with Wg an d E* determined by Eqs. ( p. 13 ), (3.43), and (1.7), and dimensionless function Fd are given by Eqs. (2.15), 
(3.52) and ( [3.53 ). We notice that the DoS singularity gets smeared by SO scattering when V > Wd- This is in contrast 
with conventional superconductivity which is known to be stable with respect to SO scattering since the latter does 
not violate T-invariance. 

This guess turns to b e corr ect. In order to demonstrate it o ne h as to show that not only the Cooperon but also 
Hikami boxes from Sec. Ill B| are modified according the rule (5.5). Taking into account the spin orbit scattering 
Eq. (5.1) in impurity lines on diagrams Figs. ^ and j^, we obtain 



B^ 2 >(u; 1 ,u; 2 ,Q 1 ,Q 2 ) 
2 



(2^ ) 3 



-i(wi + uj 2 ) + D{Q\ + Ql) + 
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+Lo 2 + cu 3 ) + D(Ql + Qj+ Ql) + 2t so ] . 
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With the same rigor as in Sec. IIIB, we conjecture that 

C„ 



B^{ Uj ,Qj} 



(2™ ) 2 



j'=i 



— J2 - iuj i + °Q. 



3Tso / 



therefore, rule ( |5.5| ) is satisfied which gives Eq. (p.7|). 

Finally we emphasize that, in Eq. ( |5.7|) complex argument of the fun ction Fd{z) shou ld be located on the physical 
sheet: z = \z\ expit/j and — 7r < ip < it. Thus, the functions \fz in Eqs. (2.15) and ( 3.5 2| ) are defined to map complex 
plane — 7r < arg(z) < it to — § < arg(z) < 5- 



B. Paramagnetic impurities, orbital effects of the magnetic field 



The derivation in the previous section suggests that any physical mechanisms of violation of either T-invariance or 
conservation of spin will have similar effect on the D oS singularity. Indeed, in OD and ID cases in the presence of 
magnetic field H and paramagnetic spins, Eq. (3.10) takes the following forma 



C(u,Q) = 



1 



i ' 



-tw + DQ 2 + 



where Ttot is a combination of phase and spin relaxation effects 



1 1 

Ttot tip 



1 



Both SO scattering and spin exchange with paramagnetic impurities (r s ) lead to spin relaxation: 



1 

T. 



(5.8) 



(5.9) 



(5.10) 



Phase may relax either due to inelastic processes or due to magnetic field effect on orbital motion of electrons 

11 1 

tin T m Th 



(5.11) 



ip iip 

When the transverse dimension of the wire (size of the grain) a, exceeds man free path /, th can be estimated a; 

2 



1 A fl 2 H „{aH\ 
— = A—^- oc D — 



(5.12) 



Here <j) = |f ~ 2 x 10- 7 Gs • cm 2 is the superconductivity magnetic flux, Et is the "transverse" Thouless energy, 
Et = D/a 2 , and fl^ is the Cooperon^ cyclotron frequency" (cyclotron frequency for a particle which mass and charge 
equal to (2D) -1 and 2e respectively)!! 



n H = 



ADeH 

he 



(5.13) 



Numerical coefficient A is not universal: it depends on the geometry of the superconductor, on the direction of 
magnetic field, etc. 



Equation ( p.8| )— ( 5.13 ) are valid also for 2D films provided H is parallel to the film plane. In this case 

a 2 fl 2 H _ D(eHa) 2 
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th 



48L> 



3c 2 r 



(5.14) 



Now, we again conjecture that Eq. ( [3.28 ) for Hikami boxes is still valid with addition of l/r fot to all DQ 2 . If so, Eq. 
(5.7) for DoS is still valid, but instead of Eq. (5^6) we should substitute 



r = 



Ttot 



(5.15) 
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Rate t v 1 in Eq. ( [5.11 ) is the contribution of inelastic collisions. This contribution will be estimated in the next 
section. Evolution of the density of states with the rate a = T/Wd for different dimensions is shown in Fig. [IB]. 
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FIG. 15. Singularity in DoS for spin-down polarized elec- 
trons for (a) OD, and (b) ID and (c) 2D systems for different 
values of dimensionless rate a = T/Wd- The DoS in two 
dimensional case is plotted for conductance g — 10. 



C. Finite temperature and inelastic processes 



Our previous consideration, strictly speaking, applies only when temperature T equals to zero. Let us now discuss 
effects of finite T. 

Temperature manifests itself through distribution of electrons in energy. This distribution substantially depends 
on T only for energies of the order of T. On the other hand, for tunneling anomalies at finite bias eV ~ E*, Eq. 
(1.7), only low temperature region T <C E* (recall that the Zeeman splitting Ez as well as superconducting gap A 
are normally of the order of E*) is of interest. Indeed, thermal broadening of the Fermi steps in the leads washes out 
any singularity as soon as T exceeds the width W of this singularity. Since the widths W d , Eqs. ( pl^ ) and ( Eg ), 
are much smaller than E* , there will be already no anomaly of the tunneling current at T ~ E* . On the other hand, 
at T <C E* the equilibrium occupation number of a state with the en ergy of the order of E* is exponentially small 
and thus can be neglected. Therefore, formally determined DoS, Eq. ( 3.50| ) is temperature independent up to terms 
of the order of exp(— E*/T). However, the observable, namely tunneling conductance axi&V), depends on T through 
the Fermi distribution np{uj/T) in the leads. To evaluate the singular part SaxieV) of the tunneling conductance, 
one has to convolute 8v (w) with the derivative of the biased Fermi distribution: 
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Sa T (eV) _ f dn F {^f-) Sv^ui) + Sv^u) 
— — / dco 



du 



Substitution of Eq. (3.5C) into Eq. (5.16) yields for <tt = <j\ + 6<tt- 



ar(eV) 
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to + aE* + iT 
W d 
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(5.17) 



To complete this discussion, let us estimate contribution of inelastic collisions of electrons, l/r„ Eq. (5.11), to Y. 
Since we are dealing with rather highly excited states (u> ~ E* ^> T), relaxation rate l/r v is determined by large 
energy transfer (~ E*), and thus is temperature independent. In all interesting cases 1/t v simply coincides with the 
energy relaxation rate l/r e . The latter (in metals, for reasonably low energies to) is determined by inelastic collision 
between electrons (see Ref. ^ for review) and can be estimated as 



1 



(5.18) 



where L w = y/D/to is the length of diffusion in time to 1 and g(L) = i/QDL d 2 is conductance of the ci-dimensional 
sample with size L. For zero dimensional case the inelastic process of this type can be neglected completely since — 

compares with the level spacing only at the energies E* of the order of the Thouless energylllii. 
Using Eqs. ( 4.13 ) we can estimate dimensionless product WdJm as 
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(5.19) 



where £ = La is the coherence length. Since A *~ CI, this estimation implies that as long as localization length Li oc 
exceeds £ the width Wd is much bigger than 1/t v and inelastic collisions are irrelevant [Li oc can be estimated from 

the equation g(Li oc ) ~ 1). 

Together with Eqs. ( 2. IS ) and ( 3.43 ) for Wd and Eqs. (5.9)— ( 5.15 ) f or T Eq. ( |5.17| ) completely describes tunneling 
anomalies at eV - ±E* for W d > T,T and d = 0, 1. For d = 2, Eq. (fyj) is valid only provided that the external 
magnetic field is parallel to the plane. The orbital effect of this perpendicular component on the tunneling anomaly 
is discussed in the next section. 



D. Magnetic field perpendicular to the film 



It was already mentioned that for zero- and one-dimensional cases orbital effects of magnetic field manifest them- 
selves through addition 1/th from Eq. (5.12), to the parameter T, see Eqs. ( |^.6[ ), (5.13), and (5.15). The same is 
true in two dimensions provided the magnetic field is parallel to the film plandj. However, effect of a component of 
magnetic field perpendicular to the plane requires a separate cousidpration. 

Similarly to usual calculation of anomalous magnetoresistanceEHla, we need to derive and to solve equation for the 
Cooperon in perpendicular magnetic field rather than to take this field into account perturbatively. This equation 
is well known to be a Schrodinger equation in imaginary time for a particle with charge 2e and mass 1/(2_D) in 
the magnetic field H . For H = cigenfunctions of this equations are plane waves. When this field H is finite 
but paral lel to the plane and weak enough (magnetic length Ijj — (hc/eH) 1 / 2 exceeds the film thickness a, or 
Eq. ( pMSj ), is smaller than "perpendicular" Thouless energy Et = D/a 2 ), it can be taken into account perturbatively. 
Eigenfunctions remain to be plane waves, so that Cooperon keeps its form Eq. (5.8), but the eigenvalues are shifted 

byTffl - 

Contrarily, even weak perpendicular component of H changes eigenfunctions, and as a result the form of the Cooper 
pole is also modified. Each eigenfunction should be characterized by the number of Landau band n and by one of the 
components of momentum Q (one of the coordinates of the guiding center) rather than by bo th components of the 
momentum. Corresponding eigenvalue equals to fln{n + 1/2), where fin is given by Eq. ( 5.13| ), i.e., it is determined 
by the number of Landau band, and it is independent of the location of the center. 

As a result, in all previous calculations DQ 2 should be substituted by Q#(n+ 1/2), and instead of integrating over 
(dQj), we have to sum over n and divide the result by Airlfj. At H = integration over dQ is limited from above by 
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DQ 2 < A, see discussion after Eq. (3.44). For the same reason we sum over n from 1 till N ~ A/S1#. In order to 
evaluate 8v a (<jS), w e have to perform calculation similar to that of Sec. [II with such changes. 
Equation ( |3.23 ) for the DoS should be rewritten as 

8v a {(J) 



where <r = ±1 corresponds to | and \ respectively, and the short hand notation 

V 



uj na = 2{u> + aE*) + i In + 

is introduced. After obvious modification of Eq. ( 3.31 ), we obtain instead of Eq. ( 3.33] ) 

Z(u>) 
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where Z(w) and m(u>) can still be expressed through (3o and /?i according t o Eq . ( 3.34 ). Consequently, functions 
f3o(uj) and /3i(w) can be connected with 7oj(C> n ) by equations similar to Eq. ( ^3.32 ) : 
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Equat ion ( 5.23 ) is valid for p = 0, 1 and for the bot h spin direc tions . We substitute Eq. ( 3.33 ) into Eqs. ( |5.20 ) and 
(5.23), car ry ou t summation over n, and obtain Eqs. (3.41) and (3.43). The scale of the singularity coincides with Wi 
from Eq. ( 3.43] ) while dimensionless function Md is substituted by 
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Here, orbital effect of the magnetic field is characterized by dimensionless parameter an — H11/W2 and ip(x) is the 
digamma function 



i>{x) = J2 
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n + 1 n + x 



-C, 



and C ~ 0.577. . . is the Euler constant. If magnetic field is weak an <C 1, we can use the asymptotic expansion 
4>(x) « In a;, x 3> 1 and recover two-dimensional result M2 from Eq. ( 3.44 ). Since function Mh depends on additional 
variable an only as on a parameter we can use the solution of Sec. Ill E| to obtain density of states: 
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wher e en ergy scale W2 is given by Eq. ( |3.43|) , we included previously discussed broadening mechanisms according to 
rule (|5.5|) with the rate T given by Eq. (|5.15|). Dimensionless function Fjj(x) is given by 
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(5.26) 



where the function y{x) is the solution of the equation 



y(x) = In 4.g - \na H - ip [ - 



-ix + y(x) 



The density of states in two dimensional films for different values of the parameter an is shown on Fig. uM 
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FIG. 16. Singularity in DoS ol two dimensional films for 
spin-down polarized electrons for for different values of di- 
mensionless magnetic field an = Qh/W2- Curves are plotted 
for conductance g — 10. 



Closing this subsection, we present the asymptotic behavior of Fjj(x ) for two limiting cases. In the weak fields 
an <C 1 the magnetic field slightly perturbs the two-dimensional result (3.53): 
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where dimensionless functions F2{x) and z(x) are defined in Eqs. ( |3 . 53| ) and ( |3.54| ) respectively. 

In the opposite limit an ^S> max(l,r /W2) the depth of the singularity is controlled by solely cyclotron frequency 
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VI. EXPERIMENTS ON Ai FILMS 

The theoretical study presented in this paper was inspired by the experimental work of Wu, Williams, and Adamsi. 
These authors studied tunneling anomalies in ultra-thin (about 4nm thick) Al films, which were driven into param- 
agnetic state by parallel magnetic field H > Hn ~ 4.8T. Both zero bias anomaly and anomalies at biases close to 
Zeeman splitting Ez were observed. The authors attempted to fit the experimental results by the theory of Ref. 
developed for normal metals and superconductors at T > T c . The agreement appeared to be reasonable with one 
important exception: the positions of the satellite singularities E** was lower in energies than that predicted by 
Ref. [|: experimentally it was fitted as 

E**~E Z -E ] E = 0.l7meV. (6.1) 



In our theory eV s — E* < Ez, see Eq. (1.7), and the discrepancy is reduced even though it does not disappear. 
For instance, minimal value of E* corresponds to the phase transition point Ez = \/2A and, according to Eq. (jOj), 
equals to A(v / 2 + l)/2 ~ 0A7meV, since A ~ 0.39meT^. Experimental value eV s ~ 0.38mey is about 20% (rather 
than 33% in comparison with Ref. smaller. We do not have enough data to speculate about possible sources 
of this discrepancy. More experiments with serious quantitative analysis are needed to verify the present theory. 
Nevertheless, it may be worthwhile to briefly discuss here how the other experimental findings of Ref. ^ compare with 
the theoretical conclusions. 

In Ref. ||, authors presented and discussed tunneling conductance G(V, H) as a function of the bias voltage and 
magnetic field for two samples. Both samples were granular Al films about Anm thick. Sheet resistance were different 
Ra^ — 4.2KQ, R^' ~ 2.0KQ. For both samples dips of the tunneling conductance at V = ±V S were observed. The 
widths at half minimum(WHM) of these dips for both samples were about 0.15 ~ 0.2meV, while the depths were 
found to be quite different: \6G/G\M ~ 0.12 and \SG/G\^ ~ 0.05. 
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One can interpret these experimental results in two different ways. The first interpretation is based on the assump- 
tion that the granular structure of the films is irrelevant, and they can be approximated as homogeneous 2D objects. 
Given the superconducting gap A ~ 0.39777,6^, Zeeman splitting Ez — 0.57meV at magnetic field H = 5T, and R a , 



one can use Eq. ( pU3|) to determine W 2 : IF 2 (1) ~ OmmeV; IF 2 (2) ~ O.OWmeV. Since WHM of the anomaly should 



be compared with approximately 2H2, see Fig. |l5|c, we hav e not great but reasonable agreement, especially for the 
first sample. However the law Wi cx \n(g)/g from Eqs. ( 3.43 ) and ( 3.53| ) seems to contradict the experiment. 

The alternative interpretation is based on the approximation of weakly connected Al grains: in the first approxima- 
tion we neglect the coupling between the grains. This allows to use 0D expression for the width of the singularity VjZg., 



see Eq. ( p. 13 ). Given the electron concentration in Al n = 1.8 x 10 23 cm 3 and their Fermi energy Ep = 11.8eV 



bare DoS.-jCstimates as v ~ 2 x 10 22 (eV cm 3 ) -1 . Assuming that the grains in lateral directions have typical size 
b ~ 3077770, and that the film thickness is a ~ 4nm, we can estimate the mean level spacing 8 ~ \/{ab 2 v) ~ 0.03meV. 
Substitution of this value of 5 into Eq. ( 2.13| ) gives Wq — O.llmeV. This is in a good agreement with the experiment, 
since WHM at zero dimensions according to Fig. [l5| a should be compared with 2Wo — 0.22. 

In order to understand the substantial difference in amplitudes of the tunneling anomalies for the two samples, let 
us discuss the effect of coupling between the grains. This coupling results in a finite dwell time Tdweii which an electron 
spends in a given grain before tunneling into a neighboring one. We can determine Tdweii from _D-the constant of the 
diffusion at times bigger than T dweU using the relation D ~ b 2 / {2T dwel {). Given the sheet resistance Rn , DoS 
and the film thickness a, one can estimate D as D^ 1 ~ 0.2cm 2 /sec, ~ 0Acm 2 /sec. As a result 



OmmeV 



J2) 



OAmeV. 



Now we can expla in the difference in the depths of the anomalies in the two samples assuming that Ti/Td contributes 
to r from Eq. (5.15). One can see from Fig. [l5| a that the dip at V = Wo/2 is approximately twice as deep as the one 
at r = Wq. At the same time, WHMs in these two cases are close to each other. 

Note that Eq. ( |5~7| ) with T — h/r dwe u can be justified only for T dwea W d > h. Theoretical investigation of the 
crossover between 0D and 2D behavior in granular films goes beyond the framework of this paper, though such a 
study can be important for a quantitative discussion of experiments. 

From the perpendicular critical ficlcLiJc^ ~ 1.5T and the dimensions of a grain one can estimate 1/th and find 
that it is irrelevant for the experimcntn. The same is correct for spin-orbit scattering. Recent studies of tunneling 
through Al grainsO show that the difference of the g-factor from 2 is very small not only in average, but also for a 
given orbital as well. Both H/th and h/r so are probably smaller than O.OlmeV and much smaller than h /Tdweii ■ 
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FIG. 17. Position of the minimum in the DoS as the func- 
tion of the Zeeman splitting Ez for A ~ 0.4m«y : (a) Theo- 
retical prediction for the superconducting stateu, see Fig. [l]a; 
(b) Our theoretical prediction, Eq. (1.7) for the paramagnetic 
state; (c) Half distance between maxima in the DoS in the 
superconducting state: (d) E* — Ez law pre dicted for the 
normal metalcl; (e) Approximation of Eq. (1.7) by a straight 
line. 



Let us now return to the discussion of the dip location. Note that a dip of DoS at finite bias exist in both 
superconducting and paramagnetic states. According to idealized Fig. |l|a in superconducting side of the Clogston- 
Chandrasekhar phase transition this anomaly is located at eV^ = A + Ez/2 (line "a" in Fig. 17). However, probably 
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due to the smearing of the DoS singularities, experimentally the minimum was found in the middle between two peaks 
in the DoS, i.e. at eV ~ A (line "c" at Fig. [FT]). As it was already mentioned, the experimentally found position of the 
singularity is lower than our theoretical prediction ([O]). In fact there were no jump in E* observed at the point of the 
first order phase transition. This discrepancy may be due to the inhomogcneous broadening of the transition-different 
granulars may have slightly different A. Another possibility is illustra ted on Fig. 
where the measurements were done the theoretical dependence Eq. ( |l.7| ) (line 
by 
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;, b » 



In the interval of magnetic fields 
at Fig. [I?]) can be approximated 



E* w rE z - 0.17meV 



(6.2) 



fit flfUj), 
than its bare value, 



with numerical factor r ~ 1.15 is slightly larger than 1 (line "e" at Fig. |l7]). Comparing Eq. (3.2) with the experimental 
we see that the theory would agree with experiment very well if we assume that actual g-factor is smaller 
i.e. g L = 2/r » 1.72. 



VII. CONCLUSION 



This paper is devoted to the anomalies of the tunneling density of states of low dimensional [d — 0, 1,2) super- 
conductors in external magnetic field. We concentrated on the Clogston - Chandrasekhar (CC) phase transition, i.e. 
the destruction of the superconductivity by the magnetic field by virtue of the Zeeman splitting. As a result normal 
paramagnetic state of electrons is formed. 

The main conclusion we can draw from our study of CC state is, that despite this state being normal (mean-field 
superconducting order parameter vanishes) , it is drastically different from a usual normal metal with some attractive 
interaction. The latter state appears, e.g., in a superconductor at temperatures higher than the transition temperature 
T c . The difference becomes apparent when one studies excited states rather than those closed to the ground state. 

Superconducting fluctuations in a usual normal disordered metal were known to contribute to the zero-bias tunneling 
anomaly as well as to Zeeman anomalies at the bias eV equal to Zeeman splittingp. However, these contributions 
(effects of the interaction in Cooper channel) are similar or weaker than effects of Coulomb repulsion of electrons, 
unless the system is anomalously close to the transition, i.e. it is not in Levanyuk-Ginzburg region. This means that 
the effects of superconducting fluctuations can be taken into account perturbatively almost everywhere, (except the 
very vicinity of the transition temperature) if the system is not too dirty. The perturbative approach (expansion in 
inverse powers of the conductance g) is valid as long as all of the characteristic length scales involved into the problem 
do not exceed the localization length Li oc . 

Tunneling anomalies in CC normal state studied by us are quite different. First of all, its position eV = E* , see 



Eq. (1.7), is different from the Zeeman splitting Ez- However, what is more important, the perturbative corrections to 
the density of states v(iS) are much singular at u> close to E* than the same order in g^ 1 corrections in usual normal 
metals. Because of this, the perturbative approach fails at paramet ricall y wider energy interval \eV — E*\ < Wd 



around the singular bias than that for the normal metal. Using Eqs. ( 3.43| ), one can check that the length scale L\y 



which corresponds to Wd is much less than , provided Li oc exceeds superconducting coherence length £. Indeed, 
since the localization length can be estimated as L^ c ~ Dv^ and L^ c ~ I exp (.Dr/ 2 )) , (where I is the mean free 
path, D is the diffusion coefficient, and is d - dimensional DoS) the characteristic spatial scale corresponding to 
the singularity, Lw d , can be written as 



J W 1 



2/3 
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The fact that Wd 3> D/Lf oc makes it necessary and also possible to go beyond the perturbation theory - one has 
to sum only most diverging terms, and it is allowed to neglect usual weak localization and interaction corrections. It 
turns out to be possible to sum directly whole series of the perturbation theory and thus determine the shapes of the 
singularities at all dimensions. 



The singularities are characterized by their widths Wd given by Eqs. (2.13) and (3.43). For zero dimensional grains 
our theory predicts a hard gap in the density of states with a given spin direction, centered at u) = E* . For one 
dimensio nal w ire the shape becomes universal (independent on v^ 1 ' and D) when energy is measured in units of Wi, 
see Eq. ( |3.52| ). It means that the depth of the anomaly is universal. In the case of two dimensio nal fil m the depths 
of the anomaly is not universal and behaves as the inverse logarithm of the conductance, see Eq. |T5l| . 



2cS 



The reason for the effects of superconducting fluctuations in CC metal to be dramatically enhanced in comparison 
with the usual case is the presence the pole-like singularity in the correlation function of these fluctuations. This pole 
at a finite frequency appears due to the fact that CC transition is of the first order. Contrarily, the temperature driven 
transition from superconductor to normal metal is of the second order, and in a usual normal state the correlator 
of the superconducting fluctuations is a smooth function of the frequency, i. e. any superconducting excitation decay 
very rapidly. We believe that the strong anomalies of the excitation spectrum at finite energies is a generic feature of 
any state created as a result of a first order quantum phase transition. 
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